Influence of Magnetism on Phonons in CaFe 2 As 2 Via Inelastic X-ray Scattering 
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In the iron pnictides, the strong sensitivity of the iron magnetic moment to the arsenic position 
suggests a significant relationship between phonons and magnetism. We measured the phonon dis- 
persion of several branches in the high temperature tetragonal phase of CaFe2As2 using inelastic 
x-ray scattering on single-crystal samples. These measurements were compared to ab initio calcu- 
lations of the phonons. Spin polarized calculations imposing the antiferromagnetic order present 
in the low temperature orthorhombic phase dramatically improve agreement between theory and 
experiment. This is discussed in terms of the strong antiferromagnetic correlations that are known 
to persist in the tetragonal phase. 



Iron-pnictide materials are currently the subject of 
enormous scientific activity. The discovery of super- 
conductivity with T c 's up to 55K and the ensuing ef- 
forts to understand the possible electronic, magnetic and 
phononic mechanisms responsible have already revealed 
intriguing information. 1,2 For example, first principles 
calculations using density functional theory for LaFeAsO 
strongly suggest that electron-phonon coupling is not suf- 
ficient to explain superconductivity. 3 Other calculations 
demonstrate strong dependence of the iron magnetic mo- 
ment on the arsenic atomic position suggesting an inter- 
play between certain phonons and magnetism Ex- 
periments and calculations of the crystal structure and 
magnetism under apressure for CaFe2As2 show a dra- 
matic sensitivity of the magnetic moment to changes in 
the atomic positions. 7 The motivation for the present in- 
vestigation was to study phonons in CaFe2As2 to ascer- 
tain any anomalies related to the structure and possible 
magnetic interactions. 

We report measurements of the phonon dispersion in 
the paramagnetic high-temperature tetragonal phase of 
single-crystal CaFe2As2 using inelastic x-ray scattering. 
Several phonon branches consisting of c-axis polarized 
As modes are observed to have energies and intensities 
in poor agreement with non-spin-polarized (NSP) band 
structure calculations of the phonon dispersion and struc- 
ture factor. However, the imposition of antiferromagnetic 
ordering in the tetragonal phase by spin-polarized (SP) 
calculations brings the dispersion into better agreement 
with the experimental data. Strong antiferromagnetic 
correlations have been observed by neutron scattering 
above the AFM ordering temperature 8 and extend above 
room temperature. 9 In addition, X. F. Wang et. alJ^ 
measured a substantial linear increase in the magnetic 



susceptibility of BaFe2As2 up to 700K. Such behavior 
may indicate spin fluctuations above the AFM ordering 
temperature.— Similar behavior has also been reported 
for CaFe2As2 up to at least 300K^ Resistivity mea- 
surements also suggest spin fluctuations and spin scat- 
tering above the AFM transition temperature^ Given 
the known sensitivity of the Fe magnetic moment to the 
As z-position in the crystal, the measurements and cal- 
culations of the phonon structure factor indicate an im- 
portant coupling between local magnetic order and inter- 
atomic force constants. 

The CaFe2As2 sample used for measurements along 
(110) was prepared by a stoichiometric elemental mixture 
in Sn metal flux. Handling of all elements was carried 
out in a N2-filled glove box and 0.2 grams of Ca, Fe, and 
As and 4.0 grams of Sn metal were used. The mixture 
was loaded into an alumina tube plugged with ceramic 
wool and sealed in a silica tube under vacuum. The mix- 
ture was heated to 1050°C for 14 hours and kept there 
for 8 hours, followed by cooling to 600° C at a rate of 20 
° c /hour . The molten Sn flux was then filtered with a cen- 
trifuge. The resulting square plate crystals of CaFe2As2 
were identified by X-ray powder diffraction. All other 
measurements were performed on samples that were pre- 
pared as described elsewhere . 14 i 15 

Inelastic x-ray scattering measurements were per- 
formed at Sector-3-ID-C of the Advanced Photon Source 
at Argonne National Laboratory. The incident energy 
was set to 21.66 keV with an energy resolution of 2.379 
meV full- width- at-half-maximum. Measurements were 
performed on a single-crystal of CaFe2As2 measuring ap- 
proximately 2x2 mm with a thickness of 100/im. The 
sample was mounted in both the (H0L) and (HHL) 
tetragonal planes in a closed-cycle Helium refrigerator 
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and most measurements were performed at room tem- 
perature, although some scans were also performed at 
lower temperatures. 

Phonons were measured using constant-Q energy scans 
along the tetragonal (0,0,5 + 0, (1,0,3 + 0, (1+^,0,3), 
and (£,£,4) directions at room temperature. The energy 
transfer scale for each scan was calibrated for monochro- 
mator and analyzer temperature drifts. The scans were 
fit to several peaks using a pseudo-Voigt line profile. 
The normalized pseudo-Voigt function is given in eqUl 
where fc (x; T)and Jl (x; T) are normalized Gaussian 
and Lorentzian functions respectively. The mixing pa- 
rameter r] = 0.393, and resolution full-width-at-half- 
maximum (FWHM) T — 2.379 rneV was determined 
from fits to the elastic scattering width of Plexiglas. 

f pV = (l-v)fG(x;T)+ V f L (x;T) (1) 

FiguresQJa) andQJb) show line scans consisting of several 
phonon excitations at Q = (0, 0, 5.75) and Q = (1.5, 0, 3) 
respectively along with fits to the data. The peak posi- 
tions for these and other scans were obtained from the fits 
and used to construct the dispersion of phonon branches 
along the different scan directions. 

In order to understand the features of the phonon dis- 
persion, the experimental measurements were compared 
to ab initio calculations of the phonons. The phonon 
dispersion was calculated using Density Functional Per- 
turbation Theory (DFPT). 16 There is not yet a consensus 
on the proper lattice parameters to use for these calcu- 
lations. In nonmagnetic calculations, relaxing the lattice 
parameters results in a large contraction of the c-axis. 
In spin-polarized calculations, the lattice distorts into an 
orthorhombic structure. With these difficulties in mind, 
experimental lattice parameters were used.— In addition, 
there is controversy over the appropriate internal Arsenic 
parameter to use. We chose relaxed positions so that all 
forces were zero. For non-magnetic and spin-polarized 
calculations, the relaxed As z-position is za s = 0.3575 
and zas = 0.3679 respectively. Mittal et. al£ use the 
experimental value of za s = 0.372(1). The pseudopo- 
tentials chosen used the Perdew-Burke-Ernzerhof (PBE) 
exchange correlation functional Settings of an 8x8x8 
k-mesh and 24 eV and 400 eV energy cutoffs for the 
wavefunctions and charge density were chosen to en- 
sure that the precision of the calculated phonon disper- 
sion was better than 1.0 meV. Due to the paramagnetic 
state of CaFe2As2 at room temperature, non-magnetic 
calculations were first performed. Phonon frequencies 
were calculated on a 2x2x2 q-mesh, and then interpo- 
lated along several symmetry directions. The resulting 
phonon frequencies and eigenvectors were used to cal- 
culate the dynamical structure factor along the selected 
scan directions. The dynamical structure factor, which is 
proportional to the x-ray scattering intensity, is given in 
eq|2p^ The Debye- Waller, Wd (Q) factor was set equal to 
zero, and the scattering length is proportional to the 
atomic number of the corresponding atom, and cr J d (q) is 
the eigenvector corresponding to the normalized motion 
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Figure 1: (a) Constant-Q line scan at Q = (0,0,5.75) . (b) 
Constant-Q line scan at Q = (1.5, 0, 3.0) . Experimental data 
are given by points and pseudo-Voigt fits by the solid line. 
Dashed red and blue lines correspond to calculations of the 
phonon structure factor for spin-polarized and non-magnetic 
calculations, respectively. 

of atom d in the j th phonon branch. 

\H j (0)\ 2 

Sj (q, = ' q , ' (1 + n 3 (q)) 5 {u - u> (q)} (2) 

H « (Q) = E -7=r exp ^ Wd (Q) + lQ ' d) i Q • °* (q) i (3) 

d vM d ^ 

For comparison with experiment, the delta functions 
were broadened with a pseudo-Voigt function whose pa- 
rameters were chosen to match the experimental resolu- 
tion. FigO] shows the comparison of the calculated struc- 
ture factor to the data for two scans. For the scans shown, 
the non-magnetic calculations show poor agreement with 
the phonons observed near ~ 17 and 20 meV. FigsE^a) 
and [3Ka) show several observed phonons (white dots) 
along (00L) and (H03) compared to the non-magnetic 
calculated phonon dispersion weighted by the structure 
factors. 

The agreement between calculations and the data 
dramatically improve when spin-polarized calculations 
were performed using a 4x4x4 k-mesh in a supercell us- 
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Figure 2: Phonon dispersion weighted by the structure factors 
along (00L) (a) without and (b) with antiferromagnetic order. 
The white dots are experimental data points. 
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Figure 3: Phonon dispersion weighted by the structure factors 
along (H03) (a) without and (b) with antiferromagnetic order. 
The white dots are experimental data points. 



ing the observed "stripe" AF structure of orthorhombic 
CaFe2As2- PBE overestimates the magnetic moment, 
with a value of 2.39 /ie/Fe versus 0.80(5) /i#/Fe observed 
experimentally. 20 Similar discrepancies between the cal- 
culated and observed magnetic moment were discussed 
in phonon studies of BaFe2As2^ Software limitations 
prevented constraining the magnetic moments while per- 
forming phonon calculations. FigO] shows that the spin- 
polarized calculations not only provide better agreement 
with the phonon energies, but even the phonon intensities 
(eigenvectors) show marked improvement. Along (00L) 
in figl2Kb), the agreement of the middle two branches 
at 17 and 20 meV is dramatically improved by includ- 
ing the AF ordering which reduces the energy splitting 
of these two branches. Examination of the eigenvec- 
tors of these modes shows that they have Ai symme- 
try and consist of longitudinally polarized Ca and As 
modes. SP calculations also show marked improvement 
for phonon branches in the same energy range along the 
(100) direction which connect to the two Ai branches at 



the zone boundary. FigOJb) shows a scan at (1.5,0,3) 
where SP calculations dramatically improve agreement 
with the two A3 phonon branches consisting of c-axis 
polarized transverse Ca and As modes in the 14-20 meV 
range. Figl3^b) also shows improved agreement of the 
SP calculations of select phonon branches along (H03) to 
the measured dispersion. Measurements of other phonon 
branches along (110) and (10L) also show better agree- 
ment with the SP calculations (not shown). 

Along both (00L) and (#00), the upper branch that 
softens with the introduction of AF order is associated 
with c-axis polarized As vibrations. Such modes have 
been predicted to have a strong influence on magnetism 
due to the sensitivity of the Fe moment on the As z- 
position^ Surprisingly, the lower branch starting at 17 
meV are primarily c-axis polarized Ca modes in the NSP 
calculations. However, the spin-polarization introduces 
a strong mixing of these two branches of identical sym- 
metry. Such mixing leads to the increased frequency and 
changes in intensity between SP and NSP calculations 
shown in figsj2]and[3l 

The improved agreement of the spin-polarized calcu- 
lations in the paramagnetic phase may seem surprising, 
however resistivity and magnetic susceptibility measure- 
ments suggest spin fluctuations and spin scattering above 
the AFM transition temperature. Inelastic neutron scat- 
tering experiments indicate substantial antiferromagnetic 
correlations that persist above T/v = 172K 8 and up to at 
least 300 K£ These correlations are observed up to high 
energies ( > 60 meV) with correlation lengths up to 20 A, 
which indicates that strong AF correlations exist between 
sizable Fe moments even above T/v. Magnetism must be 
accounted for when considering the chemical binding and 
the interatomic forces. Cooling of the sample into the AF 
ordered orthorhombic phase at 140 K appears to have lit- 
tle influence on the position or linewidth of the phonons 
along (00L). 

In summary, we have measured the phonon dispersion 
along several high symmetry directions for the param- 
agnetic high temperature tetragonal phase of CaFe2As2. 
Spin-polarized first principles calculations are in better 
agreement with the experimental results than non-spin- 
polarized calculations. The effects of large theoretical 
magnetic moments and the details of the spin fluctuations 
at room temperature are issues still under investigation. 

At the end of these investigations, Mittal, et. 
al^posted data from inelastic neutron scattering experi- 
ments and non-spin-polarized band structure calculations 
for CaFe2As2. This measurement is in excellent agree- 
ment with their published dispersion curves. Considera- 
tions of magnetic effects on the phonon force constants 
and density of states were posted by I. I. Mazin and M. 
D. Johannes 6 and T. Yildirim. 4 
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